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Abstract 

Appearance of the second coherent peak in the dynamical structure factor of 
an asymmetric spin ladder is suggested. The general arguments are confirmed by 
the first order (with respect to the asymmetry) calculation for a spin ladder with 
singlet-rung ground state. Basing on this result a new interpretation is proposed for 
the inelastic neutron scattering data in the spin gap compound CuHpCl for which 
some of the corresponding interaction constants are estimated. 

1 Introduction 

Inelastic neutron scattering is an effective method for analysis of low-energy excitations 
in low-dimensional spin systems [I]. The dynamical structure factor (DSF) obtained from 
this experiment produces an essential information about the low-energy spectrum. Sharp 
peaks of the DSF line shape correspond to coherent modes while broad bands originate 
from incoherent excitation continuums. 

Theoretical study of a spin ladder DSF was developed in the papers A strong 

antiferromagnetic rung coupling corresponds to the DSF with a single coherent-mode peak 
[2] , while for a weak coupling the line shape has only an incoherent background [3] . The 
models studied in [2], [3] are symmetric under exchange of the legs because their couplings 
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along both legs are equal to each other and the same is true for the couplings along both 
diagonals. Such requirement fails for an asymmetric spin ladder. 

The compound Cu2(C 5 Hi 2 N 2 ) 2 Cl4 (CuHpCl) was first interpreted as an asymmetric 
spin ladder [I] (with non equal couplings along diagonals). However neutron scattering 
[5], [6] revealed two coherent peaks in the DSF line shape for CuHpCl. Since this type of 
behavior does not agree with the results of [2], [3] (obtained for the symmetric case!) it 
was suggested in [6] that the magnetic structure of CuHpCl is inconsistent with the spin 
ladder model. 

In this paper we show the principal difference between excitation spectrums of sym- 
metric and asymmetric spin ladders and suggest argumentation confirming the existence 
of the second coherent peak in DSF of a asymmetric spin ladder. As an example we 
calculate the DSF for a weakly asymmetric spin ladder with singlet-rung ground state 
and produce an evidence for the second coherent peak. 

2 Hamiltonian for an asymmetric spin ladder 

The general Hamiltonian of an asymmetric spin-ladder has the following form, 

jj j^symm _|_ asymm /j\ 

where H symm = T, n H n,n+x and H asymm = £ n #nJjT- The local Hamiltonian densities 
are given by the following expressions: 

Tjsymm _ rj-rung , rjleg . Tjfrust , jjcyc fcy\ 
-"71,71+1 — n n,n+l t n n,n+l ' n n,n-\-l ' n n,n-\-\- 

where 

rj-rung 7 G G 

-"71,71+1 — •J±&l,n ' °2,n> 

rrleg _ jsymm/c q . c c \ 

n n,n+\ ~ J \\ l s l,n ' °l,n+l + °2,n ' &2,n+l)i 

Tjfrust _ jsymm/Q, c , q G "\ 

^71,71+1 — J frust ' S 2,n+1 + »2,n " Z>l,n+l), 

^n,n+l = Jc((&l,n ' Sl.fj+i) (S2.f1 ' Sa.n+l) + (Sl, n • S2,n)(Si in +i • S 2 , n +l) 

— (Sl, n • S2, n +l)(S2,n - Si jTl +i)), (3) 

and 

Tjasymm jasymm/ct o C C "\ 1 jasymm 1 o o C G \ ( A\ 

n n,n+l — J \\ l S l,n ' a l,n+l — &2,n ' ^2,77+lJ + J frust l°l,7i ' S 2,n+1 — ^2,77 ' ^l.n+lj- {^) 
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Figure 1: Schematic of the magnetic structure of an asymmetric spin ladder. 



This structure is schematically represented on the Fig.l. 

It is convenient to extract from the general H symm the " singlet-rung" part H s ~ r com- 
muting with the following operator: 

n 

where Q n = |(Si jn + S2, n ) 2 - The commutativity condition 

[H s ~ r ,Q]=0, (6) 

or in equivalent form 

[K;: +1 ,Qn + Qn + l]=0, (7) 

results to the following restriction on the interaction constants for H s ~ r 

jf^L = j;r - l j r- (8) 

According to (8) the Hilbert space for H s ~ r splits on the infinite set of eigenspaces corre- 
sponding to different eigenvalues of Q [Z],[H]: 

oo 

n = J2n m , Q\ Hm =m. (9) 

m=0 

The one-dimensional subspace Ti° is generated by the single vector 

io>=ni°>»> ( io ) 
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where |0) n is the n-th rung singlet. The following restrictions: 

j s ± - r > 2jj- r , r ± - r > -r~ r , JT r + J \\~ r > i js c ~ r , (11) 



guarantee that the state (10) is the (singlet-rung) ground state for H s ~ r . The operator 
Q has a sense of the magnon number [5] associated with H s ~ r . 
Decomposition 

jjsymm _ fjs-r _|_ ^jjsymm 

will be correct only if we put some additional restrictions on the interaction constants 
of AH symm . Postulating them in the form 

^jsymm = AJ symm = ^ AJ symm = _AJ^™. (13) 

we guarantee the uniqueness of the decomposition (12). Moreover under (13) and (4) the 
local exchange relations between AH symm , H asvmm and Q have the following form: 

{&K V ™Qn + Qn + l} = 2AH^, (14) 

r tt asymm , r) \ _ orrasymm /ir\ 

l-^Vn+l )VnTVn+lJ — ^n^+l j \±<J) 

(where { , } means anti commutator). 

As it follows from (14) the term AH symm does not mix even and odd components 
in (9). Therefore the Hilbert space H splits on two invariant subspaces, of the operator 

jjsymm 

oo oo 
-y^even _|_ -y^odd -y^even \^ s 'y_^2m q_^odd ^ q_^2m+l (16) 

m=0 m=0 

From (15) follows that H as y mm mixes Ji even and H odd , however on the sector H° its 
action is trivial. Really according to (15) iZ^^^O^ 10)^+1 have to lie in the sector with 
Qn + Qn+i = 3 that is impossible because the operator Q n has only eigenvalues and 1. 
So we have 

More detailed analysis of the 16 x 16 matrix H asymm (which represent the action of H^T™ 
on the product of ra-th and n + 1-rungs) shows that it has only three (degenerative) 

i n j i it jasymm\n T asymm\o mi r r 11 T asymm j T asymm 

eigenvalues: 0, and ±./{J^ ) z + J j rust ) ■ lheretore tor small J,, and J ^ rust 
the state |0) remains to be the ground state for H s ~ r + }j as y mm _ 

Now we may suggest the following interpretation for the appearance of the second 
coherent mode in the DSF line shape of an asymmetric spin ladder. It is known [7|- 
[9] that in the strong rung-coupling regime an excitation spectrum of a symmetric spin 
ladder has coherent modes of two types, the one-magnon triplet state lying in 7i ocM and 
three bound two-magnon states (with total spin 0,1,2) lying in '}i even _ The ground state 
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also lies in / }{ even . In the Born approximation a scattering neutron creates a new state 
by flipping a single elementary spin. It is a principal fact that the excited state lies in 
l~L odd . By this reason in the symmetric case only the subspace 7i ei,en excites during the 
the scattering process. However even a little asymmetry results to excitations from 'H even 
and in particular the bound two-magnon mode with total spin 1 which is respective for 
the appearance of the second coherent peak in the DSF. 

In the next sections we shall confirm our arguments by studying the simplest model 
for which AH symrn = and the ground state exactly has the form (10). 



3 One and two-magnon states for H s r 

The eigenstates of H s ~ r in the sectors with Q = 1 and Q = 2 may be obtained exactly 
[7J,|8J. From now we shall concern on this special model omitting the upper indexes "s— r" 
or " symm" in notation of interaction constant J c , Jii and Jf rus t- In other words we shall 
study the model (2)-(3) with additional restrictions (8) and (11) on J±, J\\, Jf rus t and J c . 
According to the following formulas: 

^- r +ll°)n|l)n + l = (\j±~ jJ C )|0) n |l>" +1 + ^|l)^|0) n+1 , 



Hn7n+l £ a/3-y\l)n\l)li+l = ( J ± ~ J \\ ~ ^c/4)£ Q /3 7 | 1)J 1)Z+1, ( 19 ) 



^-+li 1 )ni0)n + l = (-J ± --J c )|l)^|0) n+1 + ^|0) n |l)« +1 , (1* 



where a, /3, 7 = 1, 2, 3 and |1)", is the triplet associated with a n-th rung: 

|i>« = (s? iB - sy|o>, (s? >B + sy|i>2 = ie aPl \i)i (20) 

one- and (spin-1) two-magnon states for H s ~ r have the following form [Tj,[S]: 
\k,magn)% = -±= £ e lkn . . .\0) n ^\l) a n \0) , 



ft 

CO CO 

\k h k 2 ,scatt)% = Z^ tt (k 1 ,k 2 ) Yl £*^a scatt (m,n-k 1 ,k 2 )...\l)i...\l)Z..., 

m=— 00 n=m+l 

00 00 

\k,bound)% = Z&JJz) E e aPl a boun \m,n-,k)...\lf m ...\l)l.., (21) 

m=— oo n=m+l 

where 



a scatt (m, n; kt, k 2 ) = c 12 e^ kim+k2n) - c 21 e i(k2m+kin) , 

a bound (m, n; k) = e iu{m+n)+v{m - n) , u = - + (1 - rA) J. (22) 

2 AJ 2 
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Here C ab = cos - A ie *^ , A x = 5/4 - Jy/J c , u > and 

k 

cos - = |A 1 |e -1 '. 
2 1 1 

The normalization factors, 



Z scatt (ki,h) = v / 2iV,/cos2^±^-2A 1 cos^±^cos^-^ + A? ) 



N cos 2 | 
A\ - cos 2 | ' 



depend on iV the number of rungs. 

The corresponding dispersion laws are the following: 

E ma ° n (k) = J± - ^J c + J c cosk, 
E scatt (k u k 2 ) = 2J ± -3J C + J c (cosfci + cosfc 2 ), 

E bound {k) = 2J±+(Ai _ 3)Jc+ ^L cos 2^ 

As it follows from (25) the one-magnon gap E^f n and the one-magnon zone 
AE magn are given by the following formulas: 

EZ 9n = J ^~ tU - \Jcl AE ma9n = 2\J C \. 
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4 First order DSF for H s ~ r + H^ mm 

From (4) and (20) follows that, 

Trasymm h\/3h\7 _ ■( jasymm jasymm\\^\a\n,\ 

n n,n+l £ Q/37l i /nl i /n+l ~ J frust ) 1 1 /n I u /n+l 

+ i(jf— + J / 3 m )|0)n|l)^ +1 , 

SO, 

^/iV(cos 2 § - 2Ax cos f cos ^ + A ?) 



'A 2 - cos 2 \ 

«(q,magn\H as y mm \k,boundf = -2zJ^ mm (g)-^ 6 kq 5 a p, 
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where 

J 2 ' 2 

Considering H asyram as a small perturbation we may obtain the corresponding cor- 
rections for the one- and two-magnon states. In the simplest case when the one-magnon 
mode does not intersect the two-magnon sector all the first order corrections to one- 
and two-magnon dispersions vanish. First order one-magnon contributions to the S = 1 
two-magnon states are the following: 

2^2 J as v mm (kx + fc 2 ) cos S in 

\ki, k 2 , scatt)i - 



^/7V(cos 2 _ 2A1 cos cos + ^f) 
(E scatt (ki, k 2 ) -E^{k l + k 2 ))\ kl + k ^ ma ^ 



'A 2 - cos 2 \j as y mm (k) 
\k,bound) 1 = -2i A V i(jgfeomd(A;) _ Emofln(fc)) |fc,maffn) . (32) 

We use the following expression for the zero-temperature dynamical structure factor 
S aP { % u) = lim ^(OlS^l^^lS^-qJlO)^ -J5 M ), (33) 

where S(q) is the Fourier transformation of spin associated with the two dimensional 
vector q = (q,q rU ng)- Here q and q Tung are the corresponding leg and rung components. 
Since the latter has only two possible values and 7r we may study them separately, 

S(g,0) = ^e-^(S lin + S 2 , n ), S(g,vr) = ^V^S^ - S a , n ). (34) 

n n 

According to the following two formulas, 

[Q, S(q, 0)] = 0, {Q, S(q, vr)} = S(q, tt), (35) 

we may reduce the matrix elements in (33) 

(A*|S(g,0)|0> = 0, (/i|S(g,7r)|0> = ]T (n\v)(v\S(q, tt)|0>, (36) 

ueH 1 

so S a p(q,0,ui) = 0. For calculation of S a p(q,Tr,uj) let us notice that from (20) and (21) 
follows that 

(0\S a (q,n)\k,magnf = VN6 a p5 kq} (37) 
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so, the DSF has purely diagonal form, S a p(q,Tr,u>) = S a pS(q,Tr,u), where 

S(q,n,u) = ^2\(/J,\q,magn) 3 \ 2 5(uj - E^). (38) 

The unperturbed DSF corresponding only to H s ~ r consists on a single one-magnon 
coherent peak 

S {0 \q,7i,u) = 5{u-E magn {q)). (39) 

In the first order with respect to the asymmetry we have to take into account only 
the two-magnon contributions. Using the substitution 2tt — ► N f_ dk we obtain the 
following formula: 

SM{q,7T,uj) = A bound (q)^ - E bound (q)) + A scatt (q, u), (40) 

where 

4|J^— (g)| 2 (A 2 -cos 2 f) 
Abound\q) Aj(E ma 9 n (q) — E bound (q)) 2 ' 1 

4| J^ mm (g)| 2 (cos 2 1 - x 2 M)9(l - x 2 H) 
A scatt {q,uj) 7r(cos 2|_2A 1 xH + Af)(^- J E^n( g ))2- ^> 

Here Q(x) is the step function and x(u) — (uj — 2J± + 3J C )/(2J C ). 

The first term in (40) corresponds to the second coherent peak carried from the sector 

The formulas (41) and (42) will be correct only when the energy of the one-magnon 
mode is smaller than the energy of the bound state and the lower bound of the two- 
magnon continuum. Contrary due to the asymmetry mixing between Ti 1 and 7i 2 any 
intersection of the one- and two-magnon scattering sectors will result to magnon decay 
[10]. In order to avoid this possibility we shall obtain the "non-intersection" condition. 

According to (25) and (26), 

2J ± - 3J C - 2\J C \ cos h±h < E scatt (ki, k 2 ) < 2J L - 3J C + 2\J C \ cos kl ± - , (43) 

and the condition E ma9n (ki+k 2 ) < E scatt (ki, k 2 ) reduces to the following form, J c (cosk/2+ 
\J C \/(2J C )) 2 < J±/2. This inequality will be automatically satisfied for J c < 0, while for 
J c > it results to 2J± > 9J C or using (28) to an equivalent form, 

E^ n > AE magn . (44) 

The last formula has a clear interpretation. Really E™ gn — AE ma9n measures the 
difference between the one- and two- magnon sectors. When it is satisfied these sectors 
do not intersect a magnon decay is impossible and the formula (42) is correct. 
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5 Overview of experimental data for CuHpCl 



As it was suggested in [5j the compound CuHpCl corresponds to the case j® symm — o, 
^frult" 1 = J f rust- In other words it may be described by the Hamiltonian 

rrCuHpCl _ Tjrung . rrleg . rrdiag . rrcyc ( *r\ 

where the terms #™_fi, H^ 9 n+l and H^ +1 are given by (3) and 

Hn,n+1 = Jdiag&l,n&2,n+1- (46) 

XJoyo T — T _i_ jasymm _ 9 t 

iicic Udiag " frust < frust frust- 

If one suggests that the state (10) is the exact ground state then the condition (8) 
reduces to 

2J diag = J|, - ~ J c . (47) 

Under this condition (however not proved experimentally!) it is possible to estimate 
the parameters J± and J c from the formulas (25), (28) and experimental data. As it 
was presented in [TT] E™ 971 w 10.8 K, corresponds to k = it, however as it was shown in 
P2] by k = ESR measurements E gap + 2AE ma9n = 13.1 K. This data agrees with the 
neutron scattering experiments [5], [6]. From (25) and (28) follows that J c ~ 1.2 K and 
J ± ss 13.8 K. 

Unfortunately any quantitative interpretation fails for the neutron scattering data 
obtained in [6]. Really all the scans presented here correspond to the scattering with 
Qmng = 0. However as it was shown in the previous section S a p{q, 0, uS) = 0. Therefore 
appearance of the scattering peaks in Figs 9 and 10(a) of the Ref. 6 may be explained 
only by a deviation of the initial state of the ladder from the singlet-rung vacuum (10). 
The strength of this deviation may be estimated only by comparison the data presented 
in [6] with the same one related to the scattering with q rung = tt. However the latter is 
not yet obtained. 

In despite of the quantitative disagreement at q rU ng = our argumentation qualita- 
tively confirms the appearance of the second coherent peak in the structure factor. 

6 Summary 

In this paper we have demonstrated the principal difference between the excitation spec- 
trums of symmetric and asymmetric spin ladders. For the symmetric one the Hilbert 
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space splits on two invariant subspaces Ji even and H . In this case only the sector H 
gives a nonzero contribution to the dynamical structure factor. However the picture is 
quite different for an asymmetric spin ladder. The asymmetry term mixes both the sub- 
spaces and the two-magnon bound state from 'H even contribute to the DSF resulting to 
the appearance of the second coherent peak. 

As an illustration we have obtained the first order DSF for the special model of asym- 
metric spin-ladder with exact singlet rung ground state. The suggested model was applied 
to the probably asymmetric spin-ladder compound CuHpCl for which the existence of the 
second coherent peak was observed experimentally. Despite the full agreement between 
our special model and the experimental data was not confirmed some of the interaction 
constants were estimated from the inelastic neutron scattering and ESR data. 

The authors are very grateful to S. V. Maleev for helpful discussions. 
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